Resonant magnetic mode in superconducting 2-leg ladders 
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The spin dynamics of a doped 2-leg spin ladder is investigated by numerical techniques. We show 
that a hole pair-magnon boundstate evolves at finite hole doping into a sharp magnetic excitation 
below the two-particle continuum. This is supported by a field theory argument based on a SO (6)- 
symmetric ladder. Similarities and differences with the resonant mode of the high-T c cuprates are 
discussed. 
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Spin ladders materials are built from weakly coupled 
ladder units consisting of chain sub-units of spins S=l/2 
(so-called "legs") connected via some "rung" couplings 
and exhibit numerous fascinating and intriguing proper- 
ties Q . Among them, the role of the parity of the number 
of legs is remarkable: only ladders with an even number 
of legs exhibit a spin gap (ie a finite energy scale for 
triplet excitations) as seen from different characteristic 
behaviors of the spin susceptibility observed in SrCu2 03 
and SrCu3C>5 [2j, typical 2-leg and 3-leg compounds re- 
spectively. 

Besides exotic experimental properties, spin ladders 
are of special interest for theorists as well, especially the 
simple 2-leg ladder: indeed, it is believed that its ground 
state (GS) is a close realization of the Resonating Valence 
Bond (RVB) state proposed by Anderson |3| in the con- 
text of high-Tc superconductivity. The striking differ- 
ence between "even" and "odd" ladders mentioned above 
can naturally be explained in this simple picture since 
pairing nearest neighbor (NN) spins into spin singlets is 
more easily realized on rungs with an even number of 
sites. While the generic two dimensional (2D) Mott in- 
sulator (for one electron per site) is antiferromagnetic 
(AF) and may evolve into an RVB state only at finite 
doping, the 2-leg spin ladder exhibits a finite spin cor- 
relation length. It is therefore an ideal system to inves- 
tigate doping into the RVB-like spin liquid Q leading 
to d a; 2_ J( 2-like pairing As a matter of fact, pressure 
was shown to induced superconductivity in the intrinsi- 
cally doped Sr2Cai2Cu2404i+,5 ladder material hence 
bridging the apparent gap between 2-leg ladder and layer- 
based cuprates. 

In this Letter, we investigate the spin dynamics of a 
2-leg spin ladder doped with mobile holes. In 2D su- 
perconducting cuprates, a resonant mode at an energy 
around 40 meV was observed by Inelastic Neutron Scat- 
tering (INS) HQ. 

This mode was shown to be quite 
sharp both in energy and in momentum space (centered 
around the AF wavevector) and its observation seems to 
be directly linked to the appearance of superconductiv- 
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FIG. 1: Schematic picture of the low-energy triplet spectrum 
of a spin ladder doped with 2 holes. The left and right panels 
correspond to the two transverse momenta (as written on the 
plots). The two- hole continuum (h-h) and the pair-magnon 
(P-M) scattering continuum are shown as hashed regions. The 
thick line corresponds to the pair-magnon boundstate (P-M 
BS). Note that the single magnon branch is not shown (hidden 
by the P-M states). 



ity. Anomalous spectral lineshape in photoemission ex- 
periments was also interpreted as the effect of a coupling 
of the quasiparticles to a collective mode related to 
the pairing interaction. On the other hand, in ladders, 
the spin dynamics remains largely to be explored. We 
believe that a closer inspection, both experimentally and 
theoretically, of the low energy spin excitations in such a 
basic system will provide insights into the mechanism of 
pairing mediated by spin fluctuations in the 2D high-Tc 
analogs. In this study, combining different techniques, we 
provide evidences of a new magnetic mode in the two lad- 
der system and discuss similarities and differences with 
the resonant mode of the 2D cuprates. 

We shall consider here a generic 2-leg t-J ladder, 

W = Jleg Si,a ' a + J rullg Si i • Si 2 (1) 

z,a i 

+ *leg 4a C m,a + ^rung ^ 4,a C ha+l + H.C. , 
i,a i,a 

are projected hole operators and a (=1,2) 
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labels the two legs of the ladder. Isotropic couplings, 



t and Jle K = J n 



J, will be of interest 



here. A value like J = 0.4 (t is set to 1) and dopings 
between 0.1 and 0.2 are typical of the above-mentioned 
superconducting ladder materials. In this regime, the 
ladder belongs to the Luther-Emery class with a single 
zero-energy charge mode 0, ^) and the dominant su- 
perconducting correlations exhibit power-law decay |l2l | . 

Prior to the investigation of the spin dynamics at fi- 
nite doping it is very instructive to first consider a sim- 
pler case consisting of two holes (forming a hole pair) 
immersed into an infinite (or very large) ladder. Follow- 
ing previous investigations ^lj an d using finite size ex- 
trapolations , we have tentatively sketched the triplet 
excitation spectrum in Fig. ^ Apart from the two ex- 
pected types of excitations associated to (i) the dissoci- 
ation of a hole pair into two separate holes leading to a 
two-hole continuum and to (ii) magnon-hole pair scatter- 
ing states, a boundstate of the hole pair with a magnon 
state (see for details about the boundstate wavefunc- 
tion) can also be clearly identified. How this remarkable 
feature evolves at finite doping is the crucial issue we 
would like to address hereafter. 

For that purpose, the dynamic spin structure factor 
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is computed by Lanczos Exact Diagonalisations (ED) of 
2 x L periodic ladders of size L — 12 supplemented by 
a continued-fraction technique. Typical results shown in 
Fig. |21 bear similarities with the case of 2 holes in Fig.^ 
First, in the vicinity of momentum (tt,tt) a strong low 
energy incommensurate peak is observed well separated 
from the higher energy excitations. Its energy rises when 
q — * (-7T, 7r) where it becomes damped. Assuming single- 
particle excitations [TlT | at momenta ±A;f,6 [b = 1,2 for 
bonding and anti-bonding branches, kp t b ~ 7r/2), this 
feature can naturally be assigned to a fcfi + fcj?2 = 
7r(l — ri/j) spin density wave (SDW) type of excitations. 
For vanishing doping, kpi + kp 2 — ► 7r, corresponding 
to the momentum of the hole pair-magnon boundstate. 
Note also that the boundstate in Fig. ^ has a minimum 
energy of 0.34J ~ 0.14 close to the energy of the incom- 
mensurate peak in the right panel of Fig. Further- 
more, the lowest triplet excitations at q y = might be 
interpreted as the sum of two single-particle excitations 
with momenta around ±(kp,i,0) or around ±{kp£,ii) 
giving an excitation with total longitudinal momentum 
or 2kp t i ~ 2fci?2- In that case, the peaks in the left 
panel of Fig. [5] would signal the edge of this 2-particle 
continuum (or a new boundstate with small binding en- 
ergy). Note that the lowest excitation energy at q y = 
is indeed close to that of Fig. 

From the similarities between Figs. ^ and [21 ft is tempt- 
ing to interpret the q y — ir feature as a collective mode 
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FIG. 2: Dynamic spin structure factor computed on a 2 x 12 
ladder at 1/6 doping and J = 0.4. Data for and ir transverse 
momenta are shown on the left and right panels respectively, 
and q x varies from (top) to tt (bottom). The weights of the 
largest peaks (truncated and symbolized as thick arrows) are 
indicated on the plot in units where the "local" w-integrated 
weight is i ]T q S(q) = 1 - n h 
(almost) silent excitations. 
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emerging below the continuum. However, the ED data do 
not provide enough energy (and momentum) resolution 
to investigate the shape of the spectral weight in more 
details. For this purpose, the bosonized ladder model 
with SO(6)~SU(4) symmetry [3 [HI 

is of special inter- 
est. Schulz first pointed out the existence in this 
model of kink boundstates (Majorana fermions) related 
to (2kp,ir) spin density waves. Calling M the mass of 
elementary kinks 0]> the mass of these boundstates is 
M' = Ma/2 and the continuum starts at 2M. We be- 
lieve that the boundstate is, in fact, robust and more 
generic and that the above feature seen in the numer- 
ics of the t-J ladder (with lower symmetry than SO (6)) 
is a signature of it. This conjecture is partially sup- 
ported by the fact that the renormalization group flow 
of the bosonized ladder model is attracted to the direc- 
tion of SO (6) symmetry |l6|. Let us briefly summarize 
the derivation of the low energy spin dynamics within the 
SO (6) field-theoretic scheme. In bosonisation, the SDW 
operator at q x ~ 2kp (= kp t i + &f,2) and q y = n is writ- 
ten as Osdw — e l ^ p+(x ) N{x). The correlation function 
(N(x,t)N(x)) can be decomposed into[l^.l20|: 



\(0\N\B(9))\ 2 exp(M'(i- sinh 9 - r cosh 9)) 



d9 1 d9 2 \(0\N\A(9 1 )A(9 2 ))\ 2 
Mx 

exp(i (sinh 6*i +sinh# 2 ) — -Mr(cosh6>i +cosh# 2 )) 

v 

+ ... (3) 

where \B{9)) is a state which contains a single Ma- 



3 




co/M' 

FIG. 3: Contribution of the 2-kink (or 2-antikink) bound- 
states to the dynamic spin structure factor in the SO(6)- 
symmetric ladder at momentum q = (2fc.p,7r) (shown for 2 
values of K p +). The frequency is normalized to the bound- 
state energy M' — \^2M. The arrow indicates the onset of a 
small kink-kink continuum starting at ui = 2M (not shown). 
Note that the ED data of the t-J model shown on the same 
scale (setting M' as the finite size spin gap) exhibit only a 
small number of poles (5-functions). 
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FIG. 4: Spin gap of the 2-leg ladder at J = 0.4; (a) Finite 
size scaling analysis for a fixed hole number (of 2) or a fixed 
hole density computed using bosonic or boson-fermion range- 
2 CORE hamiltonians (see text) or DMRG methods. Dashed 
line are tentative fits of the data with algebraic finite size 
corrections for DMRG (~ l/£ 2 ) and exponential for CORE, 
(b) Variation of the spin gap as a function of doping. 



jorana fermion of rapidity 9, \A(6i)A(d 2 )') a state 
with 2 kinks of rapidity 9i02- Since Lorentz invari- 
ance [lH makes (0\N\B(0)) independent of 9, the contri- 
bution of the Majorana fermion to the correlation func- 
tion [HI beco mes (N(x, r)N(x)) lfeim . oxi = K (M'r/v) 
with r = x 2 + (vt) 2 . In real space, the corre- 
lation function (Osdw( x i t )Osdw(®, 0)) is given by 
{l/r) Kp+ / 2 Ko{M'r/v) where K p + is a non-universal con- 
stant whose inverse is the exponent of the power law 
decay of the superconducting correlations. The Fourier 
transform of this correlator involves a Weber Schae- 
fetlin integral[i(J and, eventually, one gets for the dy- 
namic susceptibility: xsDwiq,^) ~ Fi(l — -ftT p+ /4, 1 — 
Kp+/A; 1; (w 2 - (vq) 2 )/(M') 2 ) + where F x is the usual 
hypergeometric function and q = q x — 2k p. At the onset, 
the imaginary part of Xsdw{<1,u), i.e. S(q,u) for q = 
(2kp + q, 7r), exhibits a power-law singularity instead of 
a S- function as seen in Fig. EI In the vicinity of the onset 
one gets, XSDw(q,oj) ~ (l/(w 2 - (vq) 2 - M n )) 1 - K »+/ 2 . 
Physically, it means that, in e.g. an inelastic neutron 
scattering (INS) experiment, holons can be excited to- 
gether with the triplet boundstate, the initial neutron 
momentum being distributed among all the constituants. 
A complete characterization of such a singularity is be- 
yond the ability of present day numerical techniques but 
the results of Fig. [3 shown also in Fig. 03 revealing a low 
energy pole with a large fraction of the spectral weight 
accompanied by a few tiny peaks are consistent with such 
a scenario. 

To complete our study of the spin dynamics, we have 
combined complementary Density Matrix Renormalisa- 
tion group (DMRG) and Contractor Renormalisation 
(CORE) methods |h| supplemented by finite size anal- 



ysis to extract the collective mode energy. In the CORE 
method the ladder is decomposed into plaquette units 
and an effective interaction between adjacent plaquettes 
is constructed. In the bosonic version (B), on each pla- 
quette, only the lowest energy undoped singlet (vacuum) , 
triplet state (mapped as a triplet boson) and hole pair GS 
(mapped as a pair boson) are retained while in the more 
involved boson-fermion (BF) version the lowest one-hole 
doped plaquette states are also included. Typically, the 
effective B-hamiltonian is simply written as a sum of a 
simple bilinear kinetic term for the bosons and a quartic 
interaction. DMRG and CORE computations involving 
open and periodic boundary conditions respectively give 
rise to different types of finite size corrections as seen 
in Fig. Ufa). The resulting extrapolations of the lowest 
triplet excitation are shown in Fig. 0^b) as the function 
of the hole doping. All extrapolations with a fixed num- 
ber of 2 holes agree very well with each other and, as 
mentioned in Refs. [ill ll.'^ . the spin excitation in the 
limit of vanishing doping is not continuous due to the 
binding of the magnon with a hole pair. The main result 
shown by Fig. 0fb) is that the energy of the magnon-hole 
pair state evolves continuously, at finite doping, into the 
magnetic mode energy previously studied. Note that the 
CORE approaches (both B & BF versions) remain very 
accurate for hole doping below 10% but the bosonic ver- 
sion fails for larger doping. Physically, this signals that 
the pair size becomes too large for the spin-singlet and 
spin-triplet pair wavefunctions to be correctly captured 
by the small subset of states kept here (unless longer 
range effective interactions are included). 

Apart from singularities in S(c[,lo), the spin-triplet 
mode is expected to have effects on the spectral func- 
tions 0- We have re-examined the tunnelling density 
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FIG. 5: Tunnelling density of state of the 2-leg ladder com- 
puted on a 2 x 12 cluster. The integrated weights (sum rules) 
for electrons and holes are shown on the plot. Inset: close up 
of the low energy region. 



of states first reported in Ref. |llL looking closely at 
the low energy region. As expected, the data shown in 
Fig. El fulfill sum rules and show a quasiparticle gap Aqp 
which should correspond to the kink energy M of the 
SO (6) model. After enlarging the low energy region one 
can see two sub-peaks at an energy O ^ 0.1 from the gap 
edge. This could signal scattering of the hole or electron 
quasiparticle by a collective mode. Note that the charac- 
teristic energy f2 is indeed close to the lowest spin-triplet 
energy reported before (see Fig. 01 • 

It is interesting to compare our results to the behav- 
ior of the resonant spin excitation observed in the 2D 
copper oxide superconductors by Inelastic Neutron Scat- 
tering Remarkably, this mode was observed both 
in under-, optimally and over-doped materials and only 
below the superconducting transition temperature, Tc 
and approximately scales with Tc. While the excitation 
in the 2D cuprates occurs with wavevector (tt, 7t), an in- 
commensurate mode is found in the ladder. Although in 
2D the mode could be interpreted as a soft mode signal- 
ing the proximity of a magnetic instability , its nature 
should be different in the ladder which originates from a 
different class of "parent compound" , namely a spin liq- 
uid rather than a 2D ordered AF. Note however that, in 
2D, lower energy incommensurate spin fluctuations have 
also been seen experimentally pointing towards the exis- 
tence of a dispersive spin 1-collective mode 8j, a feature 
that might well be related to our finding in the ladder. 

To conclude, by investigating the spin dynamics of a 
doped 2-leg spin ladder, we have shown that a hole pair- 
magnon boundstate evolves into a sharp magnetic exci- 
tation at finite hole doping. Such a feature is reproduced 
by a SO (6) field theory which gives a power-law singu- 
larity at the mode energy. Analogies with the resonant 
mode of the high-T c cuprates have been discussed. 



DP & SC thank IDRIS (Orsay, France) for allocation 
of CPU-time on the NEC-SX5 supercomputer. DP also 
thanks ETH-Ziirich for hospitality and T.M. Rice for nu- 
merous discussions and comments. SRW thanks the NSF 
for support under grant DMR-0311843. SC also thanks 
A. Lauchli for help in improving some numerical codes. 



[1] E. Dagotto & T.M. Rice, Science 271, 618 (1996). 

[2] M. Azuma et al., Phys. Rev. Lett. 73, 3463 (1994). 

[3] P.W. Anderson, Science 235, 1196 (1987). 

[4] M. Sigrist, T.M. Rice and F.C. Zhang, Phys. Rev. B 49, 
12058 (1994); see also F.C. Zhang, C. Gros, T.M. Rice 
and H. Shiba, Supercond. Sci. Technol. 1, 36 (1988). 

[5] E. Dagotto, J. Riera and D.J. Scalapino, Phys. Rev. B 
45, 5744 (1992); T.M. Rice, S. Gopolan and M. Sigrist, 
Europhys. Lett. 23, 445 (1993). 

[6] M. Ueharaet al., J. Phys. Soc. J., 652 (1996); H. Mayaffre 
et al., Science 279, 345 (1998). 

[7] J. Rossat-Mignod et al., Physica (Amsterdam) 185C, 
86 (1991); H.A. Mook et al., Phys. Rev. Lett. 70, 3490 
(1993); H.F. Fong et al., Phys. Rev. Lett. 75, 316 (1995); 
P. Dai et al., Science 284^ 1344 (1999); H.F. Fong et 
al., Nature (London) 398, 588 (1999); H. He et al., 
Phys. Rev. Lett. 86, 1610 (2001). 

[8] Y. Sidis et al., cond-mat/0401328 (unpublished) and ref- 
erences therein. 

[9] Z.-X. Shen and J.R. Schrieffer, Phys. Rev. Lett. 78, 1771 
(1997). 

[10] D. Poilblanc, D.J. Scalapino, and W. Hanke, Phys. Rev. 

B, 52, 6796 (1995); C.A. Hayward and D. Poilblanc, 

Phys. Rev. B, 53, 11721 (1996). 
[11] M. Troyer, H. Tsunetsugu and T.M. Rice, Phys. Rev. B 

53, 251 (1996). 
[12] C.A. Hayward et al., Phys. Rev. Lett. 75, 926 (1995). 
[13] D. Poilblanc et al., Phys. Rev. B 62, R14633 (2000); see 

also W. Zheng et C.J. Hamer, cond-mat/0204514 
[14] S. Capponi and D. Poilblanc, Phys. Rev. B 66, 180503(R) 

(2002) and references therein. 
[15] Spectral functions for anisotropic ladders were computed 

in J. Riera, D. Poilblanc and E. Dagotto, Eur. Phys. J. B 

7, 53 (1999). 

[16] H.J. Schulz, |cond- mat/9808167 (unpublished); a 

bosonized model with a flow to a SO (6) invariant point 

is also discussed in P. Azaria, A.O. Gogolin, P. Lechem- 

inant, and A. A. Nersesyan, Phys. Rev. Lett. 83, 624 

(1999). 

[17] The total charge sector of this model decouples while 
the spin sector and the charge difference sector are de- 
scribed by the SO(6) Gross-Neveu model; D.J. Gross and 
A. Neveu, Phys. Rev. D 10, 3235 (1974). Alternatively, 
this model can be viewed as the SU(4) chiral Gross-Neveu 
model N. Andrei and J. H. Lowenstein, Phys. Rev. Lett. 
43, 1698 (1979). 

[18] These massive excitations (of mass M) carry a spin 1/2 
and are kinks of the fields characterizing (i) the total spin 
density (spinon), (ii) the difference of spin density and 
(iii) of charge current between the bonding (q y — 0) and 
antibonding (q y — tt) bands and form a 4 dimensional 
representation of SU(4). Their boundstates, the Majo- 
rana fermions, carry either spin-0 or spin-1 and form a 6 



■5 



dimensional representation of SU(4). 
[19] R. M. Konik and P. Fendley, Phys. Rev. B 66, 144416 
(2002). 

[20] F. H. L. Essler and A. M. Tsvelik, Phys. Rev. Lett. 90, 
126401 (2003). 

[21] The two-kink contribution to the correlation function for 



ui > 2M can also be computed; see E. Orignac (in prepa- 
ration) . 

[22] A.V. Chubukov, S. Sachdev and J. Ye, Phys. Rev. B 49, 
11919 (1994). 



